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I. INTRODUCTION 



The excess value of the anomalous magnetic moment of muon was reported by the E821 
collaboration at BNL [l| 

a;fP = 116 592 080(63) x 10^^^ (1) 
The Standard Model prediction for a^^^ with QED, hadronic and electroweak contributions 
is 0,3 

af^ = 116 591 785(61) x 10"^^ (2) 

with the experimental value of {g — 2)/2, the comparison gives 

Aa^, = a;fP - af^ = (295 ± 87.7) x 10~^\ (3.4a) (3) 

The 3.4 standard deviation difference between the two, may be a hint of new physics con- 
tribution. 

It has been shown that contributions from leptoquark (LQ) exchanges are capable to 
resolve the above deviation 

ABB. Leptoquarks are vector or scalar particles carrying 
both lepton and baryon numbers. LQs can be quite naturally introduced in the low-energy 



theory as a relic of a more 



scale 



unified theories (GUT) a a. In some models, it is possible to have leptoquarks at TeV 



undamental theory at some high-energy scale, such as grand 



The low-energy LQ phenomenology has received considerable attention. Possible 
LQ manifestations in various processes have been extensively investigated j^-j25|. Various 
constraints on LQ masses and couplings have been deduced from existing experimental data 
and prospects for the forthcoming experiments have been estimated. Direct searches of 
LQs as s-channel resonances in deep inelastic ep-scattering and pair production in hadron 



colliders placed lower limits on their mass Mlq > 73 — 298GeV 18| depending on the LQ 
types and couplings. The interest on leptoquarks has been renewed during the last few years 



since ongoing collider experiments have good prospects for searching these particles 27|. For 



a recent review of leptoquarks, one is referred to 



28 



Lepton Flavor Violation (LFV) are powerful tools to search for new physics. The present 
experimental limits give [isl : 

Br(/i ^ 67) < 1.2 X 10-^\ (4) 

Br(r ^ 67) < 1.1 x 10"^ (5) 

Br(r /X7) < 6.8 x 10"^ (6) 



Since effects of leptoquark interactions can manifest in a^, it is very likely that they can 



also give interesting contributions to these / l''~f processes 



6|. There are considerable 



efforts on experiments that aim at pushing the sensitivity of Br(/i — > 67) down by two order 



of magnitudes 29|. B factories and the upgraded super B factory can probe the r —>■ 67, fi'j 
decays at better sensitivities. 

The Z — >■ if decays are among the LFV interactions and the theoretica 



their branching ratios in the framework of the SM are extremely small 



30 



predictions of 



31 



results are far from the experimental limits obtained at LEPl 



3^. These 



(7) 
(8) 
(9) 

Better sensitivities are expected from the Gi ga-Z modes at future colliders, such as Interna- 
tional Linear Collider (ILC), to have 33l. l34l. l35l|: 



Br{Z - 




< 


1.7 X 10"^ 


Br{Z - 




< 


9.8 X 10"^ 


Br{Z - 




< 


1.2 X 10"^ 



Br{Z e^fi^) < 2 X 10"^ 
Br{Z e^T^) < K X 6.5 X 10" 
Br{Z fi'^T^) < K X 2.2 X 10" 



(10) 

(11) 
(12) 



with K ~ 0.2 — 1.0. It will be interesting to study the leptoquark contributions to the Z IV 
processes. 

The aim of the present paper is to study the leptoquark effects in various LFV processes 
including / — ^ /'7 and Z iV decays, while considering leptoquark contribution to as a 
solution to the muon anomalous moment discrepancy. The layout of the present paper is 
as follows: In Sec. II we introduce the formalism. We then use it in Sec. Ill to study the 
leptoquark contributions to and LFV processes including / 1''-/ and Z — > //' decays. 
Sec. IV contains our conclusions. Some formulas and low energy constraints are given in 
Appendices. 
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II. FORMALISM 



A. Scalar Leptoquark Interactions 

In this section we hst the relevant parts of the scalar leptoquark Lagrangian. We consider 
isosinglet scalar leptoquarks. The effective Lagrangian describing the leptoquark interactions 
in the mass basis is given by [lO, 24 1: 



- eQ(^^c^A,lIlY< - teQsA.Sl + leQs tan OwZ^Sl 

^^M^7^( (^3(,.-) - Q{u-)Sw)Pr - Q(u-)SwPl)uI, 

SwCw \ J 

where = 1, 2 are the indices of leptoquark, T3 = —1/2, Q^a = —2/3 are quark's isospin 
and electric charge, Qs = —1/3 is the electric charge of scalar leptoquarks Sk, a and i,j are 
quarks and leptons flavor indices and we use cw = cos^iy and sw = sm9w- The Tk,SL r 
elements of leptoquark mixing matrix that brings Sl^k to the mass basis Sk- 

= rL t'S'fc, = Tk,SuSl, (14) 



where the >S'i(/j) is the field that associates with the ejPL{R)ul terms in Clq [2J]. Note that 
in the no- mixing case (F = 1), reduce to Sl,r, which are called chiral leptoquarks, as 
they only couple to quarks and leptons in certain chirality structures. Finally, the couplings 
h and h' are 3 by 3 matrices, which give rise to various LFV processes and must be subject 
to experimental constraints. 

In this work we do not aim at a comprehensive study of the effects of all possible lepto- 
quark interactions. Instead, we try to demonstrate that a simple scalar leptoquark model 
can provide rich and interesting LFV phenomenons. 



B. Muon anomalous magnetic moment {g — 2) 



The LQ interaction is capable to generate muon anomalous magnetic moment and resolve 
the discrepancy between theoretical and experimental results. The one-loop diagrams are 
shown in Fig. [I]^a)-(b) with I = I' = fi. The extra contribution to arises from the LQ 
model due to quark and scalar leptoquark one-loop contribution is given by 
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FIG. 1: Feynman diagrams contributing to i ^ £'7 and Z — > ii' , Sk are the scalar leptoquark 
k = 1,2, u'^ are quark up with a = 1,2, 3. 



2 3 



q=l k=l 



EEx^ {\K,r,,s,\' + \K,T,,sJ'){QiunF2{x)-QsF,{x)) 



m 



K) 



m,. 



(15) 



—2/3. Our expression agrees with 
1, ...,4) depend on the variable x = 



In the above expression, Nc = 3, Qs = —1/3, Qu> 
that in [g], [l^. The kinematic loop functions Fj (z 
cslm\ are given in the appendix A. 

Using leptoquark contribution to saturate the deviation given in Eq.([3]), the leptoquark 
masses M^^ j, mixing angle O^q and couplings hql., will be constrained. 



C. £ ^ £'7 



In this subsection we give the amplitude of £ — £'7 from leptoquark exchange. According 
to the gauge invariance, the amplitude can be written as: 



ieu{p2 



7 ' 



(16) 



where is the polarization vector and q = pi — p2 is the momentum transfer. For the 
amplitude of leptoquark exchange at one-loop level as depicted in Figure. [H we have 



3 2 



^2LR 



167r2 2^ 2^ m2 

q=l k=l 



{mih'^^W,Tl^ ,Tk,SR + mi'hqehli,Tl^ ,^Tk^sL) {Q(u-)F2{x) - QsFi{x)) 



(17) 
fl8) 



with X = iTT'fu'^^/fn'sf,- The branching ratio of i ^ £'7 is: 



4r(£) 



2_R,L I I ' 



(19) 



In our numerical calculations we analyze the Brs of the decays under consideration by using 
the total decay widths of the decaying leptons T{£). 

D. Z ^ iF 



The Feynman diagrams of LFV Z decay process are shown in Fig. 1. The total contri- 
bution of all diagrams (c) and (d) can be written as 



1 



FilPr + F^j^Pl 



rrir. 



(20) 



+ 



u{pi)el{q) 



where is the Z four-momentum. The decay rates involve both F^^j^^ and F^^j^^j^j^y. 



Bt{Z 



6 Tz 
1 



T^Z 1 2 I I TpZ 1 2 



F^'rI') 



F2LRiZ)\^ + \F2RLiZ)\' 



where the form factors F^^^^ and F^^j^^p,j^^ are given by 



2LR{RL) 



(21) 



IL 



iVc 1 

167r2 



h'q,K,V^,^,Vk,SR{9sGAx) + (?,^G2(x)) (22) 
(^L - Qr) {hqih'*g'T^^^k^k,Sr.mi - /i;^/i;^,r+^ rfc,5^m^/)G'3(x) 



IR 



F,l{h^h',L^R), 



(23) 



and 



167r2 Ml 



KiKi''^Sn,k^k,SL^uA9R + 9l)G3{x) (24) 



F^RL = Fl^{h^h',L^R), 

where we have x = rnl^^/rri^g^ and the couphngs qr^l and gs are given by 



(25) 



9r 



■ {Tsiu'^) — Q{u'=) sin^ 9w) 



sin 9\y cos 9w 
9l = tan6'vy, gs = Qs^^^dw 



(26) 
(27) 



In the above expressions of Ffi(^jiy we keep only the leading term in m|/m|^. The explicit 
expressions of one loop functions G„ (n = 1, ..6) can be found in the appendix A. 



(28) 



III. NUMERICAL RESULTS AND DISCUSSION 

We are now ready to give some numerical results. The quark mass are evaluated at the 
scale of the fi = 300 GeV 26|, which is the typical leptoquark mass used in this work, 

mt = 161.4 GeV, rUc = 0.55 GeV, m„ = 11.4 x 10"^ GeV, 
and for the following quantities we use 18| 

aem = 1/137.0359, Mw = 80.45 GeV, Mz = 91.1875 GeV. (29) 
For simplicity, we assume that the couplings h and h' are real and equal to each other, i.e. 

h = h' = h*. (30) 
We use leptoquark mass splitting A = 500 GeV in our analysis, where A is defined as 



Mg^ — . Consequently, the remaining parameters in the leptoquark model are the 
mass of the light scalar leptoquark M^^, the mixing angle 6lq, and the couplings hqi. 



A. Muon Anomalous Magnetic Moment 

In this section we discuss a few phenomenological aspects of the leptoquark contributions 
to a^. In the left panel of Fig. [21 we present an scatter plot in the (M^^ — Ihqul"^) plane 



0.001 



100 200 300 400 500 600 700 800 900 1000 
Ms (GeV) 



10"' 



100 200 300 400 500 600 700 800 900 1000 
(GeV) 



FIG. 2: Scatter plot in the plane (Afg^ — l/jg^ip) in the left panel, (Mg^ —sin 26lq) in the right panel. 
These are allowed regions in the parameter space that give a^*^ = Aa^ = (295 it 87.7) x 10"^^. 

for top quark contribution (red) and charm quark contribution (green), which are allowed 
by a^Q = Aa^ = (295 ± 87.7) x 10~" [see Eq. (jS])] within the la range of data. We note 
that it is not possible to use the up quark loop contribution alone for the a^*^ = Aa^, since 
the mixing angle and couplings are strongly constrained by the vr leptonic decays (see 
Appendix B). 

In order to see the impact of the mixing angle, we present in the right panel of Fig|2] the 
allowed regions aj^'^ = Aa^ in the (M^^ — sin26'iQ) plane. We use ttem < h"^^ < 1. The 
contribution dominates around sin26'iQ ~ 0.7 both for top and charm quark contributions. 
We see that the constraint from confines the allowed range of Ms^ to Ms^ < 950 GeV 
for top quark contribution and to M^^ < 350 GeV for charm quark contribution at the 
la level. These parameter space will be used for later study of LEV processes. The light 
leptoquark mass should be below 1 TeV, if leptoquarks with couplings of electromagnetic 
strength are responsible to the deviation Aa^. It is interesting that LHC may have good 



chance to observe these particles 



27|. 



B. Lepton Flavor Violating / — > l'^ and Z W Decays 

In this section, we investigate the LEV decay processes generated by the same leptoquark 
scalar interactions. We consider only parameter space that corresponds to a^'^ = Aa^ when 
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it is appropriate. We discuss fi ^ e^y and r — > 67, yU7 decays first. 

In Fig. [3] and Fig. H] we show scatter plots of the allowed parameters in [Ms^, hqihgi') 
planes from bounds of r — fi'j, r — 67 and fj, ^ rates. Note that in the plots we use 

1.5 X 10"^^ < Br(/i 67) < 1.2 x 10-^\ 
1 X 10^^ < Br(r ^ 67) < 1.1 x 10^^ 

1 X 10"^ < Br(r ^ /i7) < 6.8 x 10"^ (31) 

where the upper bounds are from the current limits: Eqs. dl])-®, while the lower bound 
for Br(yU — > 67) is from 29| and the lower bounds for r ^ I'j are for illustration. For the 
T ^ fi'-f and fi e'-f cases the [g — 2)^ constraint is taken into account. 

For different quark contribution the couplings are bounded in the following ranges: 10~^ < 
W ^ 10-2, 10-3 < h,rKe < 1, 10-4 < htrhte < 1 and IQ-^ < hg^h.e < IQ-^. For the 
r ^ fij and /i — 67 cases, the allowed leptoquark masses are ms^ < 250 — 300 GeV 
and 1 TeV for c-quark and t— quark loop contribution, respectively. These region are 
determined from the bounds and the muon g — 2 constraint (see also Fig. [2]) at the same 
time. On the other hand the couplings governing r — > 67 decay and those generating muon 
g — 2 contribution are decoupled, the parameters corresponding to the former bounds are 
free from the latter constraint. The resulting allowed regions are larger in this cases. The 
parameters in these allowed regions will be used to predict Z IV decays. To have an idea 
of the size the allowed couplings, we give that upper bound on hqehq£r obtained form the 
present / — /'7 limits in Table [H We see that the fi e'-f constraint is more effective in 
restricting the sizes of hqihqii. 



Decay mode 




htehte' 


r ^ /i7 
r — > 67 
/X ^ 67 


< 5.29 X 10-3 

< 0.81 

< 1.45 X 10-6 


< 9.11 X 10-3 

< 0.82 

< 1.92 X 10-6 



TABLE I: Constraints on the parameters hq^hqui (q = t, c) coming from radiative FCNC processes 
induced by the scalar leptoquark using the present experimental bounds. 
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FIG. 3: Scatter plots of leptoquark parameters in {Ms^, hgghgi') planes from (£ — > £'7) bounds 
given in Eq. The left (right) figure is for the r — > /i7 (r ^ 67) case with top and charm 

quark contributions. 
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FIG. 4: Same as Fig. [3] except for the /i — > 67 case. 



In Fig. \5\ and El we give the predicted Z II' rates in correlation with Br(/ — 1'--/). 
We see that Br{Z r^e^) can reach 1.95 x 10~^, which is comparable with the present 
bound, and Br{Z jJ^T^) can reach 2.34 x 10~^, which will be accessible by future linear 
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FIG. 5: The correlation between Br(r £'7) and Br(Z ri') 'where t = e, /j,. 



10"' 



colliders. On the contrary, the current bound on the fi ^ e'y decay imposes very strong 
constraints on the related couplings as sho'wn in Table I. Hence the predicted Br{Z f^^e^) 
is rather small and is too lo^w to be observed in the near future. In Fig. [5l "we see that the 
Z — > //' rates are roughly positively correlating -with the / — ^ /'7 rates and the top quark 
loop contributions are larger than the charm quark's ones. To have observable Z r^fi^ 
and Z — i> r^e^, the r — > /i7, 67 rates are predicted to be close to the present bounds. 



In this work the analysis has been performed for the scalar leptoquark case. It is possible 
that vector leptoquarks may also contribute to {g — 2)^ and LFV processes. As shown 



in Ref. 



20|, quite often [g — 2)^ and LFV processes provide more stringent constraints 



on vector leptoquark couplings and masses than on scalar leptoquark ones. For example, 
using the measured mt and the formula given in j4|, the present Aa^ leads to a very large 
mass scale A ~ 500 TeV in the vector leptoquark case, where A was defined from the 
relation: An / = giQ/fn'j^Q. The mass scale is much larger than the corresponding mass 
scale exhibited in Fig. 2, which is found to be A ^ few - (9(10) TeV. Similarly, in / l''y 
processes, the constraints on vector leptoquark parameters are usually more severe 201]. 
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FIG. 6: The correlation between Br(/i — > 67) and Br(Z — > //e). 



IV. CONCLUSION 



Motivated by the reported discrepancy of the muon g — 2 results, we studied the lepton 
flavor violating i —>■ £'7 and Z ££' decays in the LQ model. We showed that the g — 2 
anomaly favors LQ masses in rather low-energy regime, e.g. < 1 TeV, which is within the 
reach of the forthcoming Large Hadron Collider. 

We found that leptoquarks can generate sizable LFV / — ^ ^'7 decays. The present exper- 
imental limits are used to confined the leptoquark parameter space. On the other hand, it 
is interesting to search for these LFV effects in experiments, such as MEG, B factories and 
the super B factory. 

We predict Br{Z r^e^) reaching 10~^ and Br[Z f^'^T'^) reaching 2 x 10^^, which 
can be accessible by present experiments and future linear colliders, such as ILC. On the 
contrary, the current bounds on LFV impose very strong constraints on the Br{Z fJ'^e^) 
and the ratio is too low to be observed in the near future. In this case, it is much useful to 
search for the LFV effects in — 67 decay. 
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APPENDIX A: ONE LOOP FUNCTIONS 

The loop functions Fi and Gi used in Sec. II are given by 

\2 + Zx — Qx^ + + Qx\og{x)\ 

Fiix) = ; -, (Al) 

^ ' 12(1 ' ^ ^ 



and 



Tl — 6x + Sx"^ + 2x^ — log(x j 

= ' ^^2) 

Fsix) = ^^^^ [3 - 4x + + 2 log(x)] , (A3) 

F4{x) = ^^^^^^3 [1 - + 2x \og{x)] , ( A4) 



r - 2 + - ^g^4 ^ ^^^6 _ ^2x6 log(x)l ^ ^ 

= 36(x2 - 1)^ ' 

- 36(^ 



X 



[16 - 45x2 + 36x^ - 7x^ + 12(-2 + 3x2) iog(a,)] 



3-4x2 + x^ + 41og(x) 

^3(^) = , (AT) 

„ , . 2 + 9x2 - 6x^ + x^ + ^2a;2 log(x) . . 

- 12(x2 - 1)4 ' 

^ , , 1 - 6x2 + 3x^ + 2a;6 _ i2x^ log(x) , , , 

G5(x) = ^, (A9) 

^'^''^ = 2(a:2-i)3 i-l+x'-^x' log(x)] . (AlO) 
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APPENDIX B: CONSTRAINT FORM vr ^ eue AND vr ^ iiu^ DECAYS 



0,^ 



We follow [111, |20| to constrain leptoquark parameters using pion decay data. Form the 
interactions given in Eq. f[T51) . we obtain the effective four- Fermi interaction 

h' ^'*F+ V 

il'ai"'bj^ R,k^ k,R , \ / J t-> c\ /tii \ 

Ceff = {e-PLUa){dbPRiy]) (Bl) 

•S'fe 

h h'*V^ V 

li'aiii'bj'- R,k^ k,L , \/Jr> c\ 

Mf (^*^«"") (dhPRU^) 

By using the Fierz transformation, we can rewite the Eq. OBip as 

Sk 

+ l^J^KiJ^Cj^^R,k^kAdL,bUR,a){yL,jeR,i) (B2) 

On the other hand, the conventional interation for the tt Ivi decay in the SM is given by 

£e// = -^^HM(l-75)/][rf7'^(l-75)H+h-C 
here \Vud\ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements between the con- 



stituent of the pion meson, Gi? is the 



electronic and muonic decay modes is 



36| 



^ermi couplings constant. The ratio Rth of the 



J^th 



1 + 5 



2 \ / 2 2 \ 2 

Tfit \ I ml — m ^ 



rn^, I \ ml — m'j, 

4 



= (1.2352 ±0.0001) X 10"^ (B3) 

where 5 is the radiative corrections. Thus the ratio R+h is very sensitive to non standard 
model effects (such as multi-Higges, non-chiral leptoquarks). The experimental ratio is |18l | 

Rexp = (1.2302 ± 0.004) x 10"^ (B4) 

The interference between the standard model and LQ model can be expressed by 



m 



2 



J^SM-LQ — ^th + ^th ; -7= ^ — TTo T^TTTT — T72 ] Rk'- kA^^ 

mu + md\^/2GFVudMl^me GpVudM^^ m^J ^•'^ 

At 2cr level, we get 

fm^ Rejhueh':,) m^ Re{K^,h':^) \ 
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where, 

Rmin = -1.06 X 10"^GeV"^ (B7) 
Rmax = 2.45 X lO^^GeV-^ (B8) 

The total contribution to Rsai-lq must be smaller than the differences between SM and 
experiment within the error limits allow. 
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